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Human dynamical social networks encode information and are highly adaptive. To characterize 
the information encoded in the fast dynamics of social interactions, here we introduce the entropy of 
dynamical social networks. By analysing a large dataset of phone-call interactions we show evidence 
that the dynamical social network has an entropy that depends on the time of the day in a typical 
week-day. Moreover we show evidence for adaptability of human social behavior showing data on 
duration of phone-call interactions that significantly deviates from the statistics of duration of face- 
to-face interactions. This adaptability of behavior corresponds to a different information content of 
the dynamics of social human interactions. We quantify this information by the use of the entropy 
of dynamical networks on realistic models of social interactions. Q 



INTRODUCTION 

Networks [IHS] encode information in the topology of 
their interactions. This is the main reason why networks 
are ubiquitous in complexity theory and constitute the 
underlying structures of social, technological and biolog- 
ical systems. The information encoded in social networks 
[51 [7] is essential to build strong collaborations [8] that 
enhance the performance of a society, to build reputation 
trust and to navigate [5] efficiently the networks. For 
these reasons social networks are small world [10] with 
short average distance between the nodes but large clus- 
tering coefficient. Therefore to understand how social 
network evolve, adapt and respond to external stimuli, 
we need to develop a new information theory of complex 
social networks. 

Recently, attention has been addressed to entropy mea- 
sures applied to email correspondence [llj . static net- 
works [T2TrT4] and mobility patterns [TSJ. New network 
entropy measures quantify the information encoded in 
heterogenous static networks [H1[T3]. Information theory 
tools set the limit of predictability of human mobility [15]. 
Still we lack methods to assess the information encoded 
in the dynamical social interaction networks. 

Social networks are characterized by complex organi- 
zational structures revealed by network community and 
degree correlations [IB]. These structures are sometimes 
correlated with annotated features of the nodes or of the 
links such as age, gender, and other annotated features 
of the links such as shared interests, family ties or com- 
mon work locations (TTJ [IS]. In a recent work Q35] it 
has been shown by studying social, technological and bi- 
ological networks that the network entropy measures can 
assess how significant are the annotated features for the 
network structure. 

Moreover social networks evolve on many different 
time-scales and relevant information is encoded in their 
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dynamics. In fact social networks are highly adaptive. 
Indeed social ties can appear or disappear depending on 
the dynamical process occurring on the networks such as 
epidemic spreading or opinion dynamics. Several models 
for adaptive social evolution have been proposed showing 
phase transitions in different universality classes [20lf23] . 
Social ties have in addition to that a microscopic struc- 
ture constituted by fast social interactions of the duration 
of a phone call or of a face-to-face interaction. Dynamical 
social networks characterize the social interaction at this 
fast time scale. For these dynamical networks new net- 
work measures are starting to be defined [53] and recent 
works focus on the implication that the network dynam- 
ics has on percolation, epidemic spreading and opinion 
dynamics [251129) . 

Thanks to the availability of new extensive data on 
a wide variety of human dynamics [30 34 , human mo- 
bility [13 [3BJ and dynamical social networks [37] , it 
has been recently recognized that many human activi- 
ties [5BJ are bursty and not Poissonian. New data on 
social dynamical networks start to be collected with new 
technologies such as of Radio frequency Identification De- 
vices [281 EH] and Bluetooth [31] . These technologies are 
able to record the duration of social interactions and re- 
port evidence for a bursty nature of social interaction 
characterized by a fat tail distribution of the duration 
of face-to face interactions. This bursty behavior of so- 
cial networks [TJ [5SJ I58H4T] is coexisting with modula- 
tions coming from periodic daily (circadian rhythms) or 
weakly patterns [33J- The fact that this bursty behavior 
is observed also in social interaction of simple animals, 
in the motion of rodents [33] , or in the use of words [35J , 
suggests that the underlying origin of this behavior is dic- 
tated by the biological and neurological processes under- 
lying the dynamics of the social interaction. To our opin- 
ion this problem remains open: How much can humans 
intentionally change the statistics of social interactions 
and the level of information encoded in the dynamics of 
their social networks, when they are interfacing with a 
new technology? 

In this paper we try to address this question by study- 
ing the dynamics of interactions through phone calls and 
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comparing it with face-to-face interactions. We show that 
the entropy of dynamical networks is able to quantify the 
information encoded in the dynamics of phone-call inter- 
actions during a typical week-day. Moreover we show 
evidence that human social behavior is highly adaptive 
and that the duration of face-to-face interaction in a con- 
ference follows a different distribution than duration of 
phone-calls. We therefore have evidence of an intentional 
capability of humans to change statistically their behav- 
ior when interfacing with the technology of mobile phone 
communication. Finally we develop a model in order to 
quantify how much the entropy of dynamical networks 
changes if we allow modifications in the distribution of 
duration of the interactions. 
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FIG. 1. The dynamical social networks are composed 
by different dynamically changing groups of interact- 
ing agents. In panel (A) we allow only for groups of size 
one or two as it typically happens in mobile phone commu- 
nication. In panel (B) we allow for groups of any size as in 
face-to-face interactions. 



RESULTS 
Entropy of dynamical social networks 

In this section we introduce the entropy of dynamical 
social networks as a measure of information encoded in 
their dynamics. Since we are interested in the dynamics 
of contacts we assume to have a quenched social network 
G of friendships, collaborations or acquaintances formed 
by N agents and we allow a dynamics of social interac- 
tions on this network. If two agents z, j are linked in the 
network they can meet and interact at each given time 
giving rise to the dynamical social network under study 
in this paper. If a set of agents of size N is connected 
through the social network G the agents ii, i%, . . . i n can 
interact in a group of size n. Therefore at any given 
time the static network G will be partitioned in con- 
nected components or groups of interacting agents as 
shown in Fig [T] In order to indicate that a social in- 
teraction is occurring at time t in the group of agents 
ii,Z2,...,i n and that these agents are not interacting 
with other agents, we write 9i 1 ,i 2 ,...,i n {t) = 1 otherwise 
we put Qn^.^init) — 0. Therefore each agent is inter- 
acting with one group of size n > 1 or non interacting 
(interacting with a group of size n — 1). Therefore at 
any given time 

E 0M 2 ,...,i n (*) = l. (!) 

where we indicate with Q an arbitrary connected sub- 
graph of G. The history St of the dynamical social net- 
work is given by S t = {ftj.ia, ...,«„ (t') Vi' < t}. If we in- 
dicated by p(gi 1 ,i 2 ,...,i n (t) = ljcSt) the probability that 

gi t .i 2 i n (t) = 1 given the story St, the likelihood that 

at time t the dynamical networks has a group configura- 
tion gi u i 2! ...,i„(t) is given by 

C = ]]p(g n , t2 ,..., ln (t) = l\S t y^ (2) 

Q 

The entropy S characterizes the logarithm of the typi- 
cal number of different group configurations that can be 



expected in the dynamical network model at time t and 
is given by S = — (log C) \ Sf that we can explicitly express 
as 

S = -E P ( fl ti,ta,...,i«(*) = 1 l 5 t) 1 °gP(5ti,i 2 ,...,tn(*) = M S t)- 
Q 

(3) 

According to the information theory results if the 
entropy is vanishing, i.e. S = the network dynamics is 
regular and perfectly predictable, if the entropy is larger 
the number of future possible configurations is growing 
and the system is less predictable. If we model face-to- 
face interactions we have to allow the possible formation 
of groups of any size, on the contrary, if we model the 
mobile phone communication, we need to allow only for 
pairwise interactions. Therefore, if we define the adja- 
cency matrix of the network G as the matrix a^- , the log 
likelihood takes the very simple expression given by 

£=n^w= w i(t) n p(9ii(t)=i\s t yvw 

i u|ay=l 

(4) 

with 

9i(t)+^2<lij9ij(t) = l, (5) 

3 

for every time t. The entropy is then given by 

s = -X>G*(t) = WkgpfoW = i\s t ) 

i 

-^a l]P {g l0 {t) = l|St)logpOy(i) = l\S t ). (6) 

ij 

Social dynamics and entropy of phone call 
interactions 

We have analyzed the call sequence of subscribers of 
a major euroepan mobile service provider. We consid- 
ered calls between users who at least once called each 
other during the examined 6 months period in order to 
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examine calls only reflecting trusted social interactions. 
The resulted event list consists of 633,986,311 calls be- 
tween 6, 243, 322 users. For the entropy calculation we 
selected 562, 337 users who executed at least one call per 
a day during a week period. First of all we have studied 
how the entropy of this dynamical network is affected by 
circadian rhythms. We assign to each agent i = 1,2 a 
number rii = 1,2 indicating the size of the group where 
he/she belongs. If an agent i has coordination number 
rii = 1 he/she is isolated, and if rii = 2 he/she is interact- 
ing with a group of n = 2 agents. We also assign to each 
agent i the variable t% indicating the last time at which 
the coordination number rij has changed. If we neglect 
the feature of the nodes, the most simple transition prob- 
abilities that includes for some memory effects present in 
the data, is given by a probability p n = p n (r,t) for an 
agent in state n at time t to change his/her state given 
that he has been in his/her current state for a duration 

T = t - ti. 
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FIG. 2. Mean-field evaluation of the entropy of the 
dynamical social networks of phone calls communi- 
cation in a typical week-day. In the nights the social 
dynamical network is more predictable. 

We have estimated the probability p n (T,t) in a typical 
week-day. Using the data on the probabilities p rL {r, t) we 
have calculated the entropy, estimated by a mean-field 
evaluation (Check Text SI) of the dynamical network as 
a function of time in a typical week-day. The entropy of 
the dynamical social network is reported in Fig. [2j It 
significantly changes during the day describing the fact 
that the predictability of the phone-call networks change 
as a function of time. In fact, as if the entropy of the 
dynamical network is smaller and the network is an a 
more predictable state. 



Adaptive dynamics face-to face interactions and 
phone call durations 

In this section we report evidence of adaptive human 
behavior by showing that the duration of phone calls, a 



binary social interactions mediated by technology, show 
different statistical features respect to face-to-face inter- 
actions. The distributions of the times describing human 
activities are typically broad [H H3 (HI 131 [Ml W\ , and 
are closer to power-laws, which lack a characteristic time 
scale, than to exponentials. In particular in [38] there is 
reported data on Radio Frequency Identification devices, 
with temporal resolution of 20s, showing that both distri- 
bution duration of face-to-face contacts and inter-contact 
periods is fat tailed during conference venues. 
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FIG. 3. Probability distribution of duration of phone- 
calls. (A) Probability distribution of duration of phone-calls 
between two given persons connected by a link of weight w. 
The data depend on the typical scale t*(w) of duration of the 
phone-call. (B) Probability distribution of duration of phone 
calls for people of different age. (C) Probability distribution 
of duration of phone-calls for people of different gender. The 
distributions shown in the panel (B) and (C) do not signifi- 
cantly depend on the attributes of the nodes. 



TABLE I. Typical times t*(w) used in the data collapse 
of Fig. [3j 



Weight of the link 


Typical time t*(w) in seconds (s) 


(0-2%) w max 


111.6 


(2-4%) Wmax 


237.8 


(4-8%) Wmax 


334.4 


(8-16%) w max 


492.0 


(16-32%) Wmax 


718.8 



Here we analysed the above defined mobile-call event 
sequence performing the measurements on all the users 
for the entire 6 months time period. The distribution of 
phone-call durations strongly deviates from a fat-tail dis- 
tribution. In Fig. [3]we report this distributions and show 
that these distributions depend on the strength w of the 
interactions (total duration of contacts in the observed 
period) but do not depend on the age, gender or type of 
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contract in a significant way. The distribution P M (A£j n ) 
of duration of contacts within agents with strenght w is 
well fitted by a Weibull distribution 



'(w)P w (AU 
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T*(w) 
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(7) 

with /3 = 0.47... The typical times t*(w) used for the 
data collapse of Figure 3 are listed in Table [TJ The origin 
of this significant change in behavior of humans inter- 
actions could be due to the consideration of the cost of 
the interactions (although we are not in the position to 
draw these conclusions (See Fig. [4] in which we compare 
distribution of duration of calls for people with different 
type of contract) or might depend on the different nature 
of the communication. The duration of a phone call is 
quite short and is not affected significantly by the cir- 
cadian rhythms of the population. On the contrary the 
duration of no-interaction periods is strongly affected by 
periodic daily of weekly rhythms. The distribution of 
no-interaction periods can be fitted by a double power- 
law but also a single Weibull distribution can give a first 
approximation to describe P(At no ). In Fig. [5] we report 
the distribution of duration of no-interaction periods in 
the day periods between 7AM and 2AM next day. The 
typical times r*(k) used in Figure 5 are listed in Table 
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FIG. 4. Probability distribution of duration of phone- 
calls for people with different types of contract. No sig- 
nificant change is observed that modifies the functional form 
of the distribution. 



DISCUSSION 

The entropy of a realistic model of cell-phone 
interactions 

The data on face-to-face and mobile-phone interactions 
show that a reinforcement dynamics is taking place dur- 



ing the human social interaction. Disregarding for the 
moment the effects of circadian rhythms and weakly pat- 
terns, a possible explanation of such results is given by 
mechanisms in which the decisions of the agents to form 
or leave a group are driven by memory effects dictated 
by reinforcement dynamics, that can be summarized in 
the following statements: i) the longer an agent is in- 
teracting in a group the smaller is the probability that 
he/she will leave the group; ii) the longer an agent is iso- 
lated the smaller is the probability that he/she will form 
a new group. In particular, such reinforcement principle 
implies that the probabilities p n (r, t) that an agent with 
coordination number n changes his/her state depends on 
the time elapsed since his/her last change of state, i.e., 
Pn( T ,t) — fn( T )- To ensure the reinforcement dynamics 
any function / n (r) which is a decreasing function of its 
argument can be taken. In two recent papers [TJ [41] the 
face-to-face interactions have been realistically modelled 
with the use of the reinforcement dynamics, by choosing 



fn(r) 



(t + 1) 



(8) 



with good agreement with the data when we took b n = b 2 
for n > 2 and b\ > 0, b 2 > 0. 

In order to model the phone-call data studied in this 
paper we can always adopt the reinforcement dynam- 
ics but we need to modify the probability /„(r) by a 
parametrization with an additional parameter ft < 1. In 
order to be specific in our model of mobile-phone commu- 
nication, we consider a system that consists of N agents. 
Corresponding to the mechanism of daily cellphone com- 
munication, the agents can call each other to form a bi- 
nary interaction if they are neighbor in the social net- 
work. The social network is characterized by a given 
degree distribution p(k) and a given weight distribution 
p(w). Each agent i is characterized by the size rii = 1,2 
of the group he/she belongs to and the last time ti he/she 
has changed his/her state. Starting from random initial 
conditions, at each timestep St = 1/N we take a random 
agent. If the agent is isolated he/she will change his/her 
state with probability 



/f(r) 



(r + iy 



(9) 



with r = t — ti and 61 > 0. If he/she change his/her 
state he/she will call one of his/her neighbor in the social 
network which is still not engaged in a telephone call. A 
non-interacting neighbor agent will pick up the phone 
with probability /f (r') where r' is the time he/she has 
not been interacting. 

If, on the contrary the agent i is interacting, he/she will 
change his/her state with probability f 2 '(t\w) depending 
on the weight of the link and on the duration of the phone 
call. We will take in particular 
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hgjw) 

(r + iy 



(10) 
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FIG. 5. Distribution of non-interaction times in the phone-call data. The distribution strongly depends on circadian 
rhythms. The distribution of rescaled time depends strongly on the connectivity of each node. Nodes with higher connectivity 
k are typically non-interacting for a shorter typical time scale T*(k). 



TABLE II. Typical times T*(k) used in the data collapse 
of Fig. [5] 



Connectivity 


Typical time r*(k) in seconds (s) 


k=l 


158,594 


k=2 


118,047 


k=4 


69,741 


k=8 


39,082 


k=16 


22,824 


k=32 


13,451 



where b 2 > and g(w) is a decreasing function of 
the weight w of the link. The distributions /f (r) and 
f% (t\w) are parametrized by the parameter f3 < 1. As j3 
increases, the distribution of duration of contacts and 
duration of intercontact time become broader. These 
probabilities give rise to either Weibull distribution of 
duration of interactions (if /3 < 1) or power- law distribu- 
tion of duration of interaction f3 = 1. Indeed for /3 < 1, 
the probability P2( T ) that a conversation between two 
nodes with link weight w ends after a duration r is given 
by the Weibull distribution (See Text SI for the details 
of the derivation) 

t*(w)P?(t)kWp((t + 1)/t*(w)) (11) 

with t*(w) = [2b 2 g(w)]~ 1/{1 ~ P) . This distribution well 
capture the distribution observed in mobile phone data 
and reported in Fig. [3] (for a discussion of the valid- 
ity of the annealed approximation for predictions on a 
quenched network see the Text SI ). 

If, instead of having (3 < 1 we have (3 = 1 the proba- 
bility distribution for duration of contacts is given by a 
power-law 

P^{t) oc (t + i)-[2*>2 S («0+i]_ (12) 

This distribution is comparable with the distribution 
observed in face-to-face interaction during conference 
venues [J HI]. The adaptability of human behavior, 



evident when comparing the distribution of duration of 
phone-calls with the duration of face-to-face interactions, 
can be understood as a possibility to change the exponent 
ft regulating the duration of social interactions. 

Changes in the parameter /3 correspond to a different 
entropy of the dynamical social network. Solving analyt- 
ically this model we are able to evaluate the dynamical 
entropy as a function of f3 and b\. In Fig. [6] we report 
the entropy S of the dynamical social network a function 
of f3 and b\ in the annealed approximation and the large 
network limit. In particular we have taken a network of 
size N = 2000 with exponential degree distribution of av- 
erage degree (fc) = 6, weight distribution P(w) = Cw~ 2 
and function g(w) = 1/w and b 2 = 0.05. Our aim in 
Fig. [6] is to show only the effects on the entropy due 
to the different distributions of duration of contacts and 
non-interaction periods. Therefore we have normalized 
the entropy S with the entropy Sr of a null model of 
social interactions in which the duration of groups are 
Poisson distributed but the average time of interaction 
and non interaction time are the same as in the model of 
cell-phone communication. From Fig. [6] we observe that 
if we keep b\ constant, the ratio S/Sr is a decreasing 
function of the parameter f3 indicating that the broader 
are the distribution of probability of duration of contacts 
the higher is the information encoded in the dynamics of 
the networks. Therefore the heterogeneity in the distri- 
bution of duration of contacts and no-interaction periods 
implies higher level of information in the social network. 
The human adaptive behavior by changing the exponent 
j3 in face-to-face interactions and mobile phone commu- 
nication effectively change the entropy of the dynamical 
network. 

In conclusion, in the last ten years it has been recog- 
nized that the vast majority of complex systems can be 
described by networks of interacting units. Network the- 
ory has made tremendous progresses in this period and 
we have gained important insight into the microscopic 
properties of complex networks. Key statistical proper- 
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MATERIAL AND METHODS 




FIG. 6. Entropy S of social dynamical network 
model of pairwise communication normalized with 
the entropy Sr of a null model in which the ex- 
pected average duration of phone-calls is the same 
but the distribution of duration of phone-calls and 
non-interaction time are Poisson distributed. The net- 
work size is TV = 2000 the degree distribution of the network 
is exponential with average (k) = 6, the weight distribution 
is p(w) — Cw~ 2 and g(w) is taken to be g(w) — bijw with 
62 = 0.05. The value of S/Sr is depending on the two pa- 
rameters /3,bi. For every value of 61 the normalized entropy 
is smaller for /3 — > 1. 



ties have been found to occur universally in the networks, 
such as the small world properties and broad degree dis- 
tributions. Moreover the local structure of networks has 
been characterized by degree correlations, clustering co- 
efficient, loop structure, cliques, motifs and communities. 
The level of information present in these characteristic of 
the network can be now studied with the tools of informa- 
tion theory. An additional fundamental aspect of social 
networks is their dynamics. This dynamics encode for 
information and can be modulated by adaptive human 
behavior. In this paper we have introduced the entropy 
of social dynamical networks and we have evaluated the 
information present in dynamical data of phone-call com- 
munication. By analysing the phone-call interaction net- 
works we have shown that the entropy of the network 
depends on the circadian rhythms. Moreover we have 
shown that social networks are extremely adaptive and 
are modified by the use of technologies. The statistics of 
duration of phone-call indeed is described by a Weibull 
distribution that strongly differ from the distribution of 
face-to-face interactions in a conference. Finally we have 
evaluated how the information encoded in social dynam- 
ical networks change if we allow a parametrization of the 
duration of contacts mimicking the adaptability of hu- 
man behavior. Therefore the entropy of social dynami- 
cal networks is able to quantify how the social networks 
dynamically change during the day and how they dynam- 
ically adapt to different technologies. 



In order to describe the model of mobile phone com- 
munication, we consider a system consisting of N agents 
representing the mobile phone users. The agents are 
interacting in a social network G representing social 
ties such as friendships, collaborations or acquaintances. 
The network G is weighted with the weights indicating 
the strength of the social ties between agents. We use 
Ni(to,t)dto to denote the number of agents with degree 
k that at time t are not interacting and have not inter- 
acted with another agent since time t' G (to, to + l/^0- 
Similarly we denote by N 2 ' ' w (to,t)dto the number of 
connected agents (with degree respectively k and k' and 
weight of the link w) that at time t are interacting in 
phone call started at time t' £ (to, to + 1/VY). The mean- 
field equation for this model read, 



dN*(t ,t) 
dt 



-(l + cfc)/V 1 fc (i ,i)/i(to,i) + ^2iWfe 



dN k 2 > k > w (t Q ,t) = _ 2 N^^{f ,t)h{%,t\w)+Ni^^(f)5 t t 



dt 

where the constant c is given by 



E fc '/o<*W fc, (*o.*)/i(*o>*) 
VfidtoN*' (t ,t)h(t ,t)' 



(13) 



(14) 



In Eqs. (Tl5| the rates Tr pq (t) indicate the average number 
of agentscnanging from state p = 1, 2 to state q = 1, 2 
at time t. These rates can be also expressed in a self- 
consistent way and the full system solved for any given 
choice of f\(to,t) and $2(to,t\w) (See Text SI for details). 

The definition of the entropy of dynamical social net- 
works of a pairwise communication model, is given by Eq. 
(6). To evaluate the entropy of dynamical social network 
explicitly, we have to carry out the summations in Eq. 
(6). These sums, will in general depend on the particular 
history of the dynamical social network, but in the frame- 
work of the model we study, in the large network limit 
will be dominated by their average value. In the follow- 
ing therefore we perform these sum in the large network 
limit. The first summation in Eq. (6) denotes the av- 
erage loglikelihood of finding at time t a non-interacting 
agent given a history St- We can distinguish between 
two eventual situations occurring at time t: (i) the agent 
has been non-interacting since a time t — r, and at time 
t remains non- interacting; (it) the agent has been inter- 
acting with another agent since time t — r, and at time 
t the conversation is terminated by one of the two inter- 
acting agents. The second term in the right hand side of 
Eq. (6), denotes the average loglikelihood of finding two 
agents in a connected pair at time t given a history St- 
There are two possible situations that might occur for 
two interacting agents at time t: (Hi) these two agents 
have been non-interacting, and to time t one of them de- 
cides to form a connection with the other one; (iv) the 
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two agents have been interacting with each other since a 
time t — t, and they remain interacting at time t. Tak- 
ing into account all these possibilities we have been able 
to use the transition probability form different state and 
the number of agents in each state to evaluate the en- 
tropy of dynamical networks in the large network limit 
(For further details on the calculation see the Text SI). 
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ACKNOWLEDGMENTS 



ETHICS STATEMENT 

The dataset used in this study only involved de- 
identified information and no details about the sub- 



We thank A.-L. Barabasi for his useful comments and 
for the mobile call data used in this research. MK ac- 
knowledges the financial support from EUs 7th Frame- 
work Programs FET-Open to ICTeCollective project no. 
238597 



[I] Dorogovtsev SN, Mendes JFF (2003) Evolution of net- 
works: From biological nets to the Internet and WWW. 
Oxford Univ Press. 

[2] Newman MEJ (2003) The structure and function of com- 
plex networks. SIAM Rev 45: 157-256. 

[3] Boccaletti S, Latora V, Moreno Y, Chavez M, Hwang DU, 
et al. (2006) Complex networks: Structure and dynamics. 
Phys Rep 424: 175-308. 

[4] Caldarelli G (2007) Scale-Free Networks. Oxford Univ 
Press. 

[5] Barrat A, Barthelemy M, Vespignani A (2008). Dynamical 
processes on complex networks. 

[6] Granovetter M (1973) The strength in weak ties. Am J 
Socio 78: 1360-1380. 

[7] Wasserman S, Faust K (1994) Social Network Analysis: 
Methods and applications. Cambridge Univ Press. 

[8] Newman MEJ (2001) The structure of scientific collabo- 
ration networks. Proc Natl Acad Sci USA 98: 404-409. 

[9] Kleinberg JM (2000) Navigation in a small world. Nature 
406: 845. 

[10] Watts DJ, Strogatz SH (1998) Collective dynamics of 
small-world networks. Nature 393: 440-442. 

[II] Eckmann JP, Moses E, Sergi D (2004) Entropy of dia- 
logues creates coherent structures in e-mail traffic. Proc 
Natl Acad Sci USA 101: 14333. 

[12] Bianconi G (2008) The entropy of randomized network 
ensembles. Europhys Lett 81: 28005. 

[13] Bianconi G, Coolen ACC, Perez- Vicente CJ (2008) 
Entropies of complex networks with hierarchically con- 
strained topologies. Phys Rev E 78: 016114. 

[14] Anand K, Bianconi G (2009) Entropy measures for net- 
works: Toward an information theory of complex topolo- 
gies. Phys Rev E 80: 045102. 

[15] Song C, Qu Z, Blumm N, Barabasi AL (2010) Limits of 
predictability in human mobility. Science 327: 1018-1021. 

[16] Castellano C, Fortunato S, Loreto V (2009) Statistical 
physics of social dynamics. Rev Mod Phys 81: 591-646. 

[17] Palla G, Barabasi AL, Vicsek T (2007) Quantifying social 
group evolution. Nature 446: 664-667. 

[18] Ahn YY, Bagrow JP, Lehmann S (2010) Link communi- 
ties reveal multiscale complexity in networks. Nature 466: 
761-764. 

[19] Bianconi G, Pin P, Marsili M (2009) Assessing the rele- 
vance of node features for network structure. Proc Natl 
Acad Sci USA 106: 11433-11438. 



[20] Davidsen J, Ebel H, Bornholdt S (2002) Emergence of 
a small world from local interactions: Modeling acquain- 
tance networks. Phys Rev Lett 88: 128701. 

[21] Marsili M, Vega-Redondo F, Slanina F (2004) The rise 
and fall of a networked society: A formal model. Proc 
Natl Acad Sci USA 101: 1439-1442. 

[22] Holme P, Newman MEJ (2006) Nonequilibrium phase 
transition in the coevolution of networks and opinions. 
Phys Rev E 74: 056108. 

[23] Vazquez F, Egufluz VM, San-Miguel M (2008) Generic 
absorbing transition in coevolution dynamics. Phys Rev 
Lett 100: 108702. 

[24] Tang J, Scellato S, Musolesi M, Mascolo C, Latora V 
(2009) Small-world behavior in time- varying graphs. Phys 
Rev E 81: 055101. 

[25] Holme P (2005) Network reachability of real-world con- 
tact sequences. Phys Rev E 71: 046119. 

[26] Vazquez A, Racz B, Lukacs A, Barabasi AL (2007) Im- 
pact of non-poissonian activity patterns on spreading pro- 
cesses. Phys Rev Lett 98: 158702. 

[27] Parshani R, Dickison M, Cohen R, Stanley HE, Havlin S 
(2009) Dynamic networks and directed percolation. Euro- 
phys Lett 90: 38004. 

[28] Isella L, Stehle J, Barrat A, Cattuto C, Pinton JF, et al. 
(2011) What's in a crowd? analysis of face-to- face behav- 
ioral networks. J Theor Biol 271: 166-180. 

[29] Karsai M, Kivela M, Pan RK, Kaski K, Kertesz J, et al. 
(2011) Small but slow world: How network topology and 
burstiness slow down spreading. Phys Rev E 83: 025102. 

[30] Barabasi AL (2005) The origin of bursts and heavy tails 
in humans dynamics. Nature 435: 207-211. 

[31] Eagle N, Pentland AS (2006) Reality mining: sensing 
complex social systems. Personal Ubiquitous Comput 10: 
255-268. 

[32] Rybski D, Buldyrev SV, Havlin S, Liljeros F, Makse HA 
(2009) Scaling laws of human interaction activity. Proc 
Natl Acad Sci USA 106: 12640-12645. 

[33] Malmgren RD, Stouffer DB, Motter AR, Amaral LA 
(2008) A poissonian explanation for heavy tails in e-mail 
communication. Proc Natl Acad Sci USA 105: 18153- 
18158. 

[34] Malmgren RD, Stouffer DB, Campanharo ASLO, Nunes- 
Amaral LA (2009) On universality in human correspon- 
dence activity. Science 325: 1696-1700. 



8 



[35] Brockmann D, Hufnagel L, Geisel T (2006) The scaling 
laws of human travel. Nature 439: 462-465. 

[36] Gonzalez MC, Hidalgo AC, Barabasi AL (2008) Under- 
standing individual human mobility patterns. Nature 453: 
779-782. 

[37] Onnela JP, Saramaki J, Hyvonen J, Szabo G, Lazer D, 
et al. (2007) Structure and tie strengths in mobile com- 
munication networks. Proc Natl Acad Sci USA 104: 7332- 
7336. 

[38] Cattuto C, den Broeck WV, Barrat A, Colizza V, Pinton 
JF, et al. (2010) Dynamics of person-to-person interac- 
tions from distributed RFID sensor networks. PLoS ONE 
5: ell596. 

[39] Hui P, Chaintrcau A, Scott J, Gass R, Crowcroft J, et al. 
(2005) Pocket switched networks and human mobility in 
conference environments. In: Proceedings of the 2005 
ACM SIGCOMM workshop on Delay-tolerant network- 
ing(Philadelphia, PA), pp. 244-251. 

[40] Borgnat ASP, Fleury E, Guillaume JL, Robardet C 
(2008) Description and simulation of dynamic mobility 



networks. Comp Net 52: 2842-2858. 
[41] Barrat JSA, Bianconi G (2010) Dynamical and bursty 

interactions in social networks. Phys Rev E 81: 035101. 
[42] Zhao K, Stehle J, Bianconi G, Barrat A (2010) Social 

network dynamics of face-to-face interactions. Phys Rev 

E 83: 056109. 

[43] Jo HH, Karsai M, Kertesz J, Kaski K (2011) Circadian 
pattern and burstiness in human communication activity. 
arXiv:11010377 . 

[44] Antcncodo C, Chialvo DR (2009) Unraveling the fluctu- 
ations of animal motor activity. Chaos 19: 033123. 

[45] Altmann EG, Pierrehumbert JB, Motter AE (2009) Be- 
yond word frequency: Bursts, lulls, and scaling in the tem- 
poral distributions of words. PLoS ONE 4: e7678. 

[46] Cover T, Thomas JA (2006) Elements of Information 
Theory. Wiley-Interscience. 



9 



Entropy of dynamical networks 
Supporting Informations 

K. Zhao, M. Karsai and G. Bianconi 



I. THE PROPOSED MODEL OF CELLPHONE COMMUNICATION 

A. Dynamical social network for pairwise communication 

We consider a system consisting of N agents representing the mobile phone users. The agents are interacting in 
a social network G representing social ties such as friendships, collaborations or acquaintances. The network G is 
weighted with the weights indicating the strength of the social ties between agents. To model the mechanism of 
cellphone communication, the agents can call their neighbors in the social network G forming groups of interacting 
agents of size two. Since at any given time a call can be initiated or terminated the network is highly dynamical. 
We assign to each agent i = 1,2, . .., N a coordination number n; to indicate his/her state. If m = 1 the agent is 
non-interacting, and if rii = 2 the agent is in a mobile phone connection with another agent. The dynamical process 
of the model at each time step t can be described explicitly by the following algorithm: 

(1) An agent i is selected randomly at time t. 

(2) The subsequent action of agent i depends on his/her current state (i.e. rii): 

(i) If rii = 1, he/she will call one of his/her non-interacting neighbors j of G with probability fi(U,t) where 
ti denotes the last time at which agent i has changed his/her state. Once he/she decides to call, agent j 
will be chosen randomly in between the neighbors of i with probability proportional to f\ (tj , t) , therefore 
the coordination numbers of agent i and j are updated according to the rule rii — > 2 and rij — > 2. 

(ii) If rii = 2, he/she will terminate his/her current connection with probability f2(ti,t\wij) where Wij is the 
weight of the link between i and the neighbor j that is interacting with i. Once he/she decides to terminate 
the connection, the coordination numbers are then updated according to the rule — > 1 and rij — > 1. 

(3) Time t is updated as t — > t + 1/N (initially t = 0) and the process is iterated until t = T max . 



B. General solution to the model 



In order to solve the model analytically, we assume the quenched network G to be annealed and uncorrelated. 
Therefore we assume that at each time the network is rewired keeping the degree distribution p(k) and the weight 
distribution p(w) constant. Moreover we solve the model in the continuous time limit. Therefore we always approximate 
the sum over time-steps of size St = 1/N by integrals over time. We use Ni(to,t)dto to denote the number of agents 
with degree k that at time t are not interacting and have not interacted with another agent since time t' € (to, tQ + l/N). 

Similarly we denote by ,k ,tu (io, t)dt the number of connected agents (with degree respectively k and k' and weight 
of the link w) that at time t are interacting in phone call started at time t' € (to, to + 1/N). Consistently with the 
annealed approximation the probability that an agent with degree k is called is proportional to its degree. Therefore 
the rate equations of the model are given by 

N<f(to,t)h(t ,t) - ckN*{to,t)h(to,t) + Nn k 2l {t)5 tta 

2N^ k '^(t ,t)f 2 (t ,t\w) + N7r k i'' w (t)S tt0 (15) 
where the constant c is given by 

Y Jk >fodtoN<?'{to,t)f 1 (to,t) 
Ek'k'JodtoNFito^M^t)' 



dN*(t Q ,t) 
dt 

dN^ k '- w {t ,t) 
dt 
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In Eqs. |l5| ) the rates ^ pq {t) indicate the average number of agents changing from state p — 1, 2 to state q — 1,2 at 
time t. These rates can be also expressed in a self-consistent way as 



4(i)4E f dt f 2 (t ,t 



\w)N k - k '- w (t ,t) 



-12 (*) = ^ 

where the constant C is given by 



dt I dt^N^t^Nl'it'^hih^hit'^ik + k 1 ) 



C = J2 [ dt k'N*'{t Q ,t)h(t ,t). 
,., Jo 



(17) 



(18) 



The solution to Eqs. (151 is given by 



Nf(to,t) = NK k 2l {t a )e- (1+ck)S to h[ta > t)dt 
N k ' k '' w {t ,t) = Nn k i k '- w (to)e~ 2l t» h{ta - Mw)dt 



(19) 

which must satisfy the self-consistent constraints Eqs. ( 17) and the conservation of the number of agents with different 
degree 



dt [N k (t ,t) + N k - k '> w (t ,t)] = Np(k). 



(20) 



In the following we will denote by P k (to,t) the probability distribution that an agent with degree k is non-interacting 
for a period from to to t and by P^(to,t) the probability that a connection of weight w at time t is active since time 
to. It is immediate to see that these distributions are given by the number of individual in a state n = 1,2 multiplied 
by the probability of having a change of state, i.e. 

P k {t Gl t) = {l + ck)f 1 {t Ql t)N k {t ()l t) 

P?(t ,t) = 2fc(t ,t\w)Y,N k ' k '> w (t ,t). (21) 



k,k' 



C. Stationary solution with specific fi(to,t) and f2(to,t) 



In order to capture the behavior of the empirical data with a realistic model, we have chosen 



/ 1 (to,t) = /i(r) 



f2(t ,t\w) 



(1 

Mt\w) = 



b2g{w) 



(1 



(22) 



In Eqs. 
t-to ). 



(|22j), t is the duration 
The functions of /i(r) 



with parameters b± > 0, 62 > 0, < f3 < 1 and arbitrary positive function g(w) 
time elapsed since the agent has changed his/her state for the last time (i.e. r = 
and /2(t|w) are decreasing function of their argument t reflecting the reinforcement dynamics discussed in the main 
body of the paper. The function g(w) is generally chosen as a decreasing function of w, indicating that connected 
agents with a stronger weight of link interact typically for a longer time. We are especially interested in the stationary 
state solution of the dynamics. In this regime we have that for large times t 3> 1 the distribution of the number of 
agents is only dependent on t. Moreover the transition rates Tr pq (t) also converge to a constant independent of t in 
the stationary state. Therefore the solution of the stationary state will satisfy 

N k (t Q ,t) = N k (r) 



N, 



k,k ,w 1 



T k.k ,1 



(r) 



'(t ,t) = N« 

Kpqit) = n pq . (23) 
The necessary condition for the stationary solution to exist is that the summation of self-consistent constraints given 



by Eq (16[) and Eq, 
Eqs. (23). The convergence depends on the value of the parameters 6 



( 18 ) together with the conservation law Eq. 



( 20 ) converge under the stationary assumptions 
b\, j3 and the choice of function g{w). In 
particular, when < j3 < 1, the convergence is always satisfied. In the following subsections, we will characterize 
further the stationary state solution of this model in different limiting cases. 
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1. Case < P < 1 



The expression for the number of agent in a given state N*(t) and N 2 ' ' w (t) can be obtained by substituting Eqs 



(22) into the general solution Eqs. (19 1, using the stationary conditions Eqs. (23). In this way we get the stationary 
solution given by 

■\rk.k' ,w / \ AT k.k'w 2b 29(™) m _m _j_ r '\i-^i k,k' ,w w ( \ (c\a\ 

N 2 { T )^ N ^i2 ' e i-e 1 J =A/V 12 m 2 (r). (24) 



To complete the solution is necessary to determine the constants 7r|x and 7Ti 2 in a self-consistent type of solution. To 
find the expression of 7r 12 '™ as a function of 7T 21 we substitute Eqs. (24) in Eq.(17) and we get 

1 r Ft 

k,k ,w 



(*) = ^liPH 







fc / dt m1(to,t)fi(t ,t) I dt' Ni' {t' ,t)h{t' ,t) 



o 



+ fc'/ dt m1(t 0> t)fi(t ,t) I d«J,JVf ' (*£,, t) 



(25) 



Finally we get a closed equation for 71-^ by substituting Eq.(25) in Eq.(20) and using the definition of c and C, given 
respectively by Eq. (16) and Eq. (18). Therefore we get 



'21 



mi(r)c?T 



P(w) / m%{T)dTdw 



x(ck m k 1 (T)f 1 ( T )dT+ mf(r)/i(T)d 



p(fc). 



(26) 



Performing explicitly the last two integrals using the dynamical solution given by Eqs. (24), this equation can be 
simplified as 



n 21 



mi(r)dT 



P(w) / m™{ T )dTdw 



p(k). 



Finally the self-consistent solution of the dynamics is solved by expressing Eq. ( 16 ) by 

= E^li(i + cfc)- x 

E^a+cfc)- 1 ' 



(27) 



(28) 



Therefore we can use Eqs. (27) and (28) to compute the numerical value of 7r 21 and c. Inserting in these equations 
the expressions for /i(t), f2(r\w) given by Eqs. (22) and the solutions (t), jV 2 ' k ' w (t) given by Eqs. (24) we get 

k h(l + ck) M^ (1+T) i-g 

1 [ ' a + T y 

(29) 



Ff{j) oc 2 M^) e -W(^)^, 

2 W (1 +T)/ 3 



The probability distributions Pi(t) and P^t), can be manipulating performing a data collapse of the distributions, 
i.e. 



r 2*H 



AiXi _,9 e i-f 



1-/3 



(30) 



with r*(fc) and t 2 (u>) defined as 



^(A)= [^(l + cfc)] »- 
r|H= [2b 2 g(w)]-^ 



(31) 



where ^! and A 2 are the normalization factors. The data collapse defined by Eqs. (30) of the curves Pf (t), P 2 1 "(t) 
and are both described by Weibull distributions. 
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FIG. 7. Data collapse of the simulation of the proposed model for cell phone communication. In the panel (A) we plot the 
probability P™(t) that in the model a pair of agents with strenght w are interacting for a period r and in the panel (B) we 
plot the probability P* (r) that in the model an agents of degree k is non-interacting for a period r The simulation data on 
a quenched networks are compared with the analytical predictions (solid lines) in the annealed approximation. The collapses 
data of Pa"(r) is described by Weibull distribution in agreement with the empirical results found in the mobile phone data. 



D. Comparisons with quenched simulations 



To check the validity of our annealed approximation versus quenched simulations, we performed a computer simu- 
lation according to the dynamical process on a quenched network. In Fig. [7] we compare the results of the simulation 
with the prediction of the analytical solution. In particular in the reported simulation we have chosen j3 = 0.5, 
bi = 0.02, hi — 0.05 and g(w) = w~ 1 , the simulation is based on a number of agent N — 2000 and for a period of 
Tmax = 10 5 , finally the data are averaged over 10 realizations and the network is Poisson with average (k) = 6 and 
weight distribution p(w) cx w~ 2 . In Fig. [TJ we show evidence that the Weibull distribution and the data collapse 
of P™(t) well capture the empirical behavior observed in the mobile phone data (Fig. 3]). The distribution of the 
non-interaction periods P*(t) in the model is by construction unaffected by circadian rhythms but follow a similar 
data collapse as observed in the real data (Fig. [5]) . The simulated data are also in good agreement with the analytical 
prediction predicted in the annealed approximation for the parameter choosen in the figure. As the network becomes 
more busy and many agents are in a telephone call, the quenched simulation and the annealed prediction of Pi(t) 
differs more significantly. 



1. Case /3 = 



For /? = the functions /i(t) and /2(t|w) given by Eqs.(22) reduce to constants, therefore the process of creation 
of an interaction is a Poisson process and no reinforcement dynamics is taking place in the network. Assigning (3 = 
to Eqs. (19), we get the solution 



JVftr) = N^ 1 e- bl(1+ck)l 



N. 



k,k' ■ 



(r) = Nir^ k ^ e - 2b29(u,)T . 



(32) 



and consequently the distributions of duration of given states Eqs. (21 1 are given by 



P%(t) cx e - 2b29{w)T 



(33) 



Therefore the probability distributions Pf (r) and PftV) are exponentials as expected in a Poisson process. 
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2. Case /? = 1 



In this section, we discuss the case for (3 = 1 such that /* (r) oc (1 + r) 1 and f 2 (T\w) cx (1 + r) 1 . Using Eqs. 



(15) we get the solution 



N*{t) = ATt4(1 + r )- bl ( 1+cfc ) 



and consequently the distributions of duration of given states Eqs. (21 ) are given by 



P 2 w (r) oc 7r&*'' w (l + t)- 262 ^- 1 . 



(34) 



(35) 



The probability distributions are power-laws. This result remains valid for every value of the parameters b\,b 2l g(w) 
(See Ref. [T] for a full account of the detailed solution of this model) nevertheless the stationary condition is only 
valid for 



&i(l + dfe) > 1 
2b 2 g(w) > 1. 



(36) 



Indeed this condition ensures that the self-consistent constraits Eqs. (16), (18) and the conservation law Eq. (20) 
have a stationary solution. 



E. Solution of the mean-field model on a fully connected network 

Finally, we discuss the mean-field limit on the model in which every agent can interact with every other agent. In 
this case, social network is a fully connected network. Therefore we use Ni(to,t) and A^toi*) to denote the number 
of agents of the two different states respectively and the rate equations are then revised to 



dt 

dN 2 (t 0: t) 
dt 



-2JVi(to,t)/i(to,t) + Mr 2 i(t)atto 
-2N 2 (t ,t)f 2 (t ,t) + Nir 12 (t)6 tto 



(37) 



Since we will refer to this model only in the framework of a null model, we will only discuss the case in which the 
dynamics of the network is Poissonian, i.e. when 



fi(t ,t) = h 
f 2 (t ,t) = b 2 . 

The stationary solution of this model is given by exponentials, i.e. 

JVi(r) =NTT 21 e- 2blT 
N 2 (t) = NTT 12 e- 2b2T . 



Finally the distributions of duration of given states expressed by Eqs. (21) are given by 

Pi(t) oc e- 2blT 



P 2 {t) oc e 



-2b 2 r 



(38) 



(39) 



(40) 



which are exponential distributions as expected in a Poisson process. 



14 



II. ENTROPY OF THE DYNAMICAL SOCIAL NETWORKS 

A. Entropy of the dynamical social networks of pairwise communication 

The definition of the entropy of dynamical social networks of a pairwise communication model, is given by Eq. (6) 
of the main body of the article that we repeat here for convenience, 

S = -^P( 5i (t) = l\S t ) log P( 9l (t) = l\St) 

i 

- *ijP(9a(t) = l\St) hgP( 9ij (t) = l\S t ) (41) 

ij 

In this equation the matrix is the adjacency matrix of the social network and gij(t) — 1 indicates that at time t the 
agents i and j are interacting while <7;(i) = 1 indicates that agent i is non-interacting. Finally S t = {gi(t'), gij(t') Vt' < 
t} indicates the dynamical evolution of the social network. In this section, we will evaluate the entropy of dynamical 
social networks in the framework of the annealed model of pairwise communication explained in detail in the previous 
section of this supplementary material. To evaluate the entropy of dynamical social network explicitly, we have to 



carry out the summations in Eq. (41 1. These sums, will in general depend on the particular history of the dynamical 
social network, but in the framework of the model we study, in the large network limit will be dominated by their 
average value. In the following therefore we perform these sum in the large network limit. The first summation in 



Eq. (41) denotes the average loglikelihood of finding at time t a non- interacting agent given a history S t - We can 
distinguish between two eventual situations occurring at time t: (i) the agent has been non-interacting since a time 
t — t, and at time t remains non- interacting; (ii) the agent has been interacting with another agent since time t — r, 
and at time t the conversation is terminated by one of the two interacting agents. In order to characterize situation 
(i) we indicate by P^-i( T ) the probability that a non-interacting agent with degree k in the social network, that 
has not interacted since a time t, doesn't change state. Similarly, in order to characterize situation (ii), we indicate 

by J^' 1 " (t) the probability that a connected pair of agents (with degrees k and k' respectively, and weight of the 
link w) have interacted since time r and terminate their conversation at time t. Given the stationary solution of the 

pairwise communication model, performed in the annealed approximation, the rates P/L^t) and P^L,{ U '(t) are given 
by 



PUir) = 1- ^ - ^ E / ^'(r') fl (r')dr' 

k' J 

Mr) 



1 - (1 + ck)- 



N 



P^ir) = (42) 



where the constant C is given by 



C = J2j k'Ni\r')h{T')dr' (43) 



and /i(t) and f2(r\w) are given in Sec. IC The variable iVf (r) indicates the number of agents of connectivity k 



nonintcracting since a time r. This number can in general fluctuate but in the large network limit it converges to 



its mean-field value given by Eq. (24 1 The second term in the right hand side of Eq. (41 1, denotes the average 
loglikelihood of finding two agents in a connected pair at time t given a history S t . There are two possible situations 
that might occur for two interacting agents at time t: (Hi) these two agents have been non-interacting, and to time t 
one of them decides to form a connection with the other one; (iv) the two agents have been interacting with each other 

since a time t — r, and they remain interacting at time t. To describe the situation (Hi), we indicate by Pfj_f 2 ( T j T ') 
the probability that two non interacting agents, isolated since time t — r and t — t' respectively, interact at time t. 
In order to describe situation (iv), we denote by P^V^ 7 ") the probability that two interacting agents, in interaction 
since a time t — r, remain interacting at time t. In the framework of the stationary annealead approximation of the 
dynamical network these probabilities are given by 

P^{r,r') = f -m^l{k + k') 



P 2 ^(r) = l- 2 ^. 
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Therefore, the entropy of dynamical social networks given by Eq. (41 1 can be evaluated in the thermodynamic limit, 
and in the annealed approximation, according to the expression 



k J ° 

pOO 

- / N^'^^P^'r^logP^'r^dr 



k.k' ,w 



IE 



k.k' 



OO p oo 



N*(t)N* (r')P^ 2 (r, r') log P^ 2 (r, r')drdr' 



k,k' 



2 E / N%> k - W (T)P^ W (T)l0gP^ 2 > W (T)dT, 
k,k'.w • J ° 



(45) 



with TVf (t) and N 2 ' ' (t) given in the large network limit by Eqs. (24) 



B. Entropy of the null model 

To understand the impact of the distribution of duration of the interactions and of the distribution of non-interaction 
periods, we have compared the entropy S of the pairwise communication model with the entropy Sr of a null model. 



Here we use the exponential mean-field model described in Section IE as our null model. In this model the agents 
are embedded in a fully connected networks and the probability of changing the agent state does not include the 
reinforcement dynamics. In fact we have that the transition rates are in dependent of time (j3 = 0) and given by 
/«(r) = 6f and f 2 R (r) 



b R . 



Following the same steps used in Sec. 



II A 



for evaluating S in the model of pairwise 



communication on the networks, it can be easily proved that the entropy Sr of the dynamical null model is given by 



57 



r 26f" 




r 26f- 






log 


dr 


N _ 




N _ 





(r)— log— dr 

2 V ) N S N 



OO POO 



N«(T)N«(T>)^l0g^LdTdT> 



^0 

oc 



26£ 
N 



NC R 
log 



1 



2b§ 



N 



NC R 
dr 



(46) 



where the constant C R is given by 



C R 



N R (r)dr, 



(47) 



and where Ni, N 2 are given , in the large network limit by their mean- field value given by Eq.(39|. In order to build 
an appropriate null model for the pairwise communication model parametrized by (/3, b\,b 2 ) ,we take the parameters 
of the null model b R and b 2 such that the proportion of the total number of agents in the two states (interacting or 
non-interacting) is the same in the pairwise model of social communication and in the null model. In order to ensure 
this condition we need to satisfy the following relation 



Jo°°N R (T)di 



E 



k.k' 



\r)dr I^N R (r)d, 



(48) 



In particular we have chosen b R — b\ and we have used Eq. ( 48 ) to determine b R 
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III. MEASUREMENT OF THE ENTROPY OF A TYPICAL WEEK-DAY OF CELL-PHONE 

COMMUNICATION FROM THE DATA 

In this section we discuss the method of measuring the dynamical entropy from empirical cellphone data as a 
function of time t in a typical weekday. This analysis gave rise to the results presented in figure 2 in the main body 
of the paper. We have analyzed the call sequence of subscribers of a major European mobile service provider. We 
considered calls between users who at least once called each other during the examined 6 months period in order to 
examine calls only reflecting trusted social interactions. The resulted event list consists of 633.986.311 calls between 
6.243.322 users. For the entropy calculation we selected 562.337 users who executed at least one call per a day during 
a working week period. Since the network is very large we have assumed that the dynamical entropy can be evaluate 
in the mean-field approximation. We measured the following quantities directly from the sample: 

• N±(t, t) the number of agents in the sample that at time t are not in a conversation since time t — r; 

• N calls (T 7 t) the number of agents in the sample that are not in a conversation since time t — r and make a call 
at time t; 

• N called (T,t) the number of agents in the sample that are not in a conversation since time t — r and are called 
at time t; 

• M m (r, t) the number of agents that at time t are in a conversation of duration r with another agent in the 
sample; 

• M out (t, t) the number of agents that at time t are in a conversation of duration r with another agent outside 
the sample; 

• M end (T,t) the number of calls of duration r that end at time t. 

Using the above quantities, we estimated the probability p calls (r,t) that an agent makes a call at time t after a non- 
interaction period of duration r, the probability p called (T,t) that an agent is called at time t after a non-interaction 
period of duration r and the probability 7t(t, t) that a call of duration r ends at time t, according to the following 
relations 

K ' ' JVi(r,t) 

T\Tcalled( j-\ 
called, .\ _ jV \ T i l ) 

P iVk(r,t) 

- M end (r,t) 
n{ - T ' ' ~ M in (r, t)/2 + M mt (t, t) ' ^ ' 

Since the sample of 562.337 users we are considering is a subnetwork of the whole dataset constituted by 6.243.322 
users, in our measurement, an agent can be in one of three possible states 

• state 1: the agent is non- interacting; 

• state 2: the agent is in a conversation with another agent of the sample; 

• state 3: the agent is in a conversation with an agent outside the sample. 



Therefore , to evaluate the entropy of the data, we can modify Eq.(41 1 into 



S{t) - - yP(9i(t) = l\S t ) log P{ 9l (t) = l\S t ) 



i 

E 

ij 



aijPigijit) = l\S t )\ogP( gij (t) = l\S t ) 



-Y^P(9'M = l|5 t )logP(^(t) = l\S t ) (50) 

i 

where et^- is the adjacency matrix of the quenched social network, gi(t) — 1 indicates that the agent i is in state 1, 
9ij(t) = 1 indicates that the agent is in state 2 interacting with agent j and g[{i) = 1 indicates the agent i is in state 
3. Finally St — {gi(t'), gij(t') g[{t) Vt' < t} indicates the dynamical evolution of the social network. To explicitly 
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evaluate Eq. (50) in the large network limit where we assume that the dependence on the particular history are 
vanishing, we sum over the loglikelihood of all transitions between different states using the same strategy in Sec. 2, 
which is 

s(t) = - x>i ( t ' *) p i^i ( r ' *) lo § p i->i ( r ' *) 

T 

- J2 Min (^ t)P*+l{T, t) log P 2 ^(T, t) 

T 

- M out (r, t)P 3 ^ (r, t) log P 3 ^(r, t) 

T 

t,t' 

T 

- ]T A^r, t)P^ 3 (T, t) log P^3(r, f) 

r 

- ]T M ollt (r, t)P 3 ^ 3 (r, t) log F 3 ^3 (t, t) . 

r 

where the probabilities of transitions between different states are given by 

Pi-»i(r,i) = 1 -p CQ " s (r,t) -p called (T,t) 

P 2 ^l(T,t)=P 3 ^l(T,t)=ir(T,t) 



P^ 2 (T,T>,t) = 



(1-7) 
C 



(T, t)p Called (T / , t) + P CaUS {T', t) P CaUea {T, t) 



calls (I 



xlled i 



P2-+2(T,t) =P 3 ^ 3 {T,t) =l-7T(T,t) 



P 1 _, 3 (r,i)=7 



P CallS {T,t) +p Call£d {T, t) 



and where C is given by 



C = ^iYi(r,i) J5 ca " ed (r,t). 



(51) 



(52) 



(53) 



Finally in 52 we have introduced a parameter 76 [0,1] to denote the portion of the calls occurring between an agent 
in the sample and an agent out of the sample. For simplicity, we assume that 7 is a constant. Substituting Eq.(52) 
into Eq.(51 ), we have performed the summation over r to obtain the value of entropy as a function of t presented in 
Figure 2 of the main body of the paper where we have taken 7 = 0.8, consistently with the data. 
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